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A description of scalar charged particles, based on the Feshbach - Villars formalism, is proposed. 
^— ■») , Particles are described by an object that is a Wigner function in usual coordinates and momenta 

' and a density matrix in the charge variable. It is possible to introduce the usual Wigner function for 

a large class of dynamical variables. Such an approach explicitly contains a measuring device frame. 
From our point of view it corresponds to the Copenhagen interpretation of quantum mechanics. 
It is shown how physical properties of such particles depend on the definition of the coordinate 
operator. The evolution equation for the Wigner function of a single-charge state in the classical 
limit coincides with the Liouville equation. Localization pecuharities manifest themselves in specific 
constraints on possible initial conditions. 
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O ■ I- INTRODUCTION 

(N 

' The question regarding the nature of the wavcfunction has its origin in the early years of quantum mechanics. For 
a long time it has been considered as a philosophical rather then a physical question |p]-^. However, it is now very 
real because of the recent theoretical and experimental progress in quantum information . 

One of the significant points in understanding the nature of the wave function is the Einstein - Podolsky - Rosen 
paradox and the existence of quantum correlations related to it. Since 1980, such specific behavior of quantum 
systems has been confirmed many times in experiments It is very important that such correlations "spread" 
' in space instantly. Nevertheless, if one adheres to the Copenhagen interpretation, there is no violation of causality 
^ ' principle. 

However, due to the fact that collapse of the wavcfunction takes place at a given instant (or maybe in a small time 
^ interval) in the whole space, one can speak of violation of causality principle and a conflict between quantum mechanics 
• I— I and special relativity in some interpretations of quantum mechanics. Following the idea of Bell and Eberhard 
that a consistent description should contain a preferred frame, relativistic classical and quantum mechanics is built in 



^ in such a way that the relativity principle is generalized and the causality is not violated. A basic assumption 

- - of this theory is that transformation from one frame to another is realized by operators which are isomorphic to the 
Lorentz group and depend on a certain four-vector (as parameter) in such a way that the instant time hyper-plain 
is invariant. This vector is interpreted as a relative to an observer four-velocity of the preferred frame. Under this 
assumption the postulate about the light speed constancy is changed to the postulate about the constancy of the 
average speed of light on closed path. 

One of the remarkable peculiarities of this theory is in the existence of a well-defined position operator which 
coincides with the Newton - Wigner position operator pT| , p^ in the preferred frame. It means that in this approach 
measurement of position does not create a particle - antiparticle couple since the odd part responsible for appearance 
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of a different charge state superposition is absent. Hence, there exists a supposition that locahzation of relativistic 
particles can bear information about presence of the preferred frame in the Universe. 

The above arguments show the fundamental role that determination of the position operator structure can play. 
Moreover, the question of whether it is a one-particle, or, in principle, a many-particle operator, is of specific interest 
as well. Therefore, it is very important to find situations where the odd part of the position operator could manifest 
itself. It would be possible to realize such observation on strongly localized (near the Compton wavelength) states of 
a separate particle. However, because such states are a problem for laboratory experiments, it is of great interest how 
this peculiarity manifests itself in a many-particle system. 

The Weyl - Wigner - Moyal (WWM) formalism is a convenient method to describe both one-particle and many- 
particle quantum systems in non-relativistic theory jl^-jl^. Nevertheless, attempts to generalize it for the relativistic 
case lead to a number of problems. The first problem is that the Weyl rule does not include time as a dynamical 
variable, and the scalar product in the Hilbert space of states is formulated for functions square integrable not over 
the whole space-time but in the three-dimensional space or in a space-like hyper-surface only. In [18] this problem was 
resolved by generalization of the spatial integration over the whole space-time without Weyl rule application. The 
WWM formalism in the framework of the stochastic formulation of quantum mechanics, where the scalar product is 
formulated for square integrable functions in the whole space-time, also leads to Lorenz-invariant expressions ]l9| . 

The matrix- valued Wigner function formalism has been developed for the general case of many component equations 
and, in particular, for 1/2 spin particles on the basis of the usual Weyl rule ||2^,^. Certainly, such equations are not 
Lorenz invariant, though average values coincide with their analogues in the usual approach. 

The next essential problem in relativistic WWM formalism is in absence of a well-defined position operator. Unlike 
the works mentioned above, where the Wigner function is determined by means of the usual position operator, in 
p2| , p3| a formalism with the Newton - Wigner position operator is developed. This approach is not Lorenz invariant 
either, and respective results differ from the standard ones. However, they can be related to [Q, where the consistent 
definition of the position operator is possible. 

The goals of this paper are in formulating the WWM relativistic formalism for scalar charged particles under the 
approach of po| and finding a set of specific peculiarities in relativistic quantum system behavior, that are related to 
the non-trivial structure of the position operator. It turns out that values of some observables depend directly on the 
position operator definition, that, as mentioned above, can be a consequence of existence of the preferred frame in 
the Universe. 

In Sec.|| we introduce the Weyl rule for matrix- valued observables in the case of scalar charged particles and discuss 
peculiarities of correspondence between classical and quantum dynamical algebras. Sec. Ill is devoted to the matrix- 
valued Wigner function and quantum Liouville equation. Here we also discuss the absence of Lorenz invariance of 
this approach, and how it can be related to the Copenhagen interpretation of quantum mechanics. In Sec, IV we 
consider how the usual Weyl rule turns into the Feshbach - Villars representation. Among the whole set of dynamical 
variables, we separate a special class of observables with the Weyl symbols independent of charge variable (charge- 
invariant variables). They have several remarkable properties. In particular, we find the relationship between their 
even and odd parts. The usual (not matrix- valued) Wigner function can be introduced for such observables. It 
is considered in Sec.^ (for a brief abstract of this approach see [Q). This object includes four components: one 
corresponds to a particle, second one corresponds to an antiparticle (even part of Wigner function), and two more 
are interference terms (odd part of Wigner function) . The evolution equation for the odd part becomes zero in the 
classical limit, and the equation for the even part coincides with its analogue in the Newton - Wigner position operator 
approach |^,^. The difference reveals itself in peculiarities of constraints on initial conditions that arc considered 
in Sec.lvl 



II. WEYL RULE AND SPECIFIC PROPERTIES OF DYNAMICAL ALGEBRA FOR SCALAR 

CHARGED PARTICLES 

In this and the following Sections we apply the methods developed in to the Klein - Gordon equation that 

is written using the Feshbach - Villars formalism ||l2| . This makes it possible to take the charge variable into account 
explicitly. With account of the fact that the Hilbert space of states for scalar charged particles has an indefinite metric, 
we have to distinguish between covariant and contravariant basis vectors and coordinates for subspacc corresponding 
to the charge variable. Hence, arbitrary state is expanded on basis vectors in the following way: 

i*)= E (1) 

Q = ±l 

where 
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(2) 



Here s = ... is an arbitrary dynamical variable and the summarizing symbol can be interpreted as an integration 
with respect to s. Greek indices take values ±1. r^ap plays the role of metric tensor (see |12 ). In notations like 
, Ti^l3, T2af3, T^ap, T^a^,..., symbols 6, Ti do not corrcspoud to specific operators and mean only values of the 
corresponding matrix elements. There are some peculiarities in notations of operators as well. A is a full operator 
that acts on all dynamical variables. A is an operator that acts on the charge variables only (or, in other words, it is 
a c- numeric matrix). 

For a consistent development of the WWM formalism we shall formulate the Weyl rule bringing into correspondence 
matrix-valued classical variables (Weyl symbols) to quantum-mechanical operators. Here we should take into account 
that in our case classical variables are operators acting on the charge variable. However, because the operator of 
quasi-probability density in the coordinate representation is proportional to the identity matrix in charge space, the 
succession of this operator and a matrix-valued Weyl symbol does not matter. Choosing it in an arbitrary way, one 
can write the Weyl rule as follows: 



-t-oo 

Aj^J2 I A^P{p,q)W^^p,q)dpdq, 

, Li 



(3) 



7=±1 



where Aa^{p,q) is the matrix-valued Weyl symbol, Aa^ stands for corresponding quantum-mechanical operator, 
Wa^{p,q) is the operator of quasi-probability density that is the Fourier transform of the displacement operator. 
Following one can obtain the expansion of Wa^{p, q) on eigenvectors of the position operator: 



q + ^)Sa'^ exp 



h 



dQ{q-^ 



(4) 



where d = 3 is dimensionality of the physical space. 

Substituting (0) to (0) one can find the expression for an arbitrary operator expansion: 



Aa 



- 



1 



1 + -^) Aj{p, q)exp 



n 



I Q 
dQdpdq ( 9 ~ 



(5) 



Similarly, such expression can be written by expending the operator of quasi-probability density on the eigenvectors 
of the momentum operator: 



A 1^ - 



1 



+ OC 



P+-^ )Aa'^{P, 9)exp 



I P 

dPdpdq ( P - Y 



(6) 



Matrix elements of this operator can be found in the form: 

If 1 

(gi|^a'^|g2) - -^^^ j A„^(p, -(91+92)) exp 



^(91 - 92 )P 



dp. 



(7) 



By changing variables in this expression in a standard way and performing Fourier transformation one can obtain the 
expression that reconstructs matrix- valued Weyl symbol through operator: 



+00 



Aj{p,q)^ / Iq 



Q 



A 1^ 



dQ. 



(8) 



In ||2^ the integral form of the Klein - Gordon equation with pseudo-differential symbols in the Feshbach - Villars 
representation (where Hamiltonian matrix has a diagonal form) is obtained. In fact, it means that the Newton - 
Wigner coordinate is used there instead of the usual coordinate. Here we obtain another integral form of the Klein - 
Gordon equation using the usual position operator: 
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+ CXD 



E 

/3=±1_ 



'^p{y)dpdy, 



(9) 



where 



(T3'- 



■ IT2 



p .^ {p~eA{x)f 
2m 



(10) 



is the matrix- valued Weyl symbol of the Hamiltonian in the general case (for a particle in electromagnetic field), which 
depends on both coordinate and momentum. Unlike ||2^ this equation takes into account that the position operator 
mixes states with different charge signs. 

Until now, consideration of the matrix-valued WWM formalism has been very similar to the usual one. Considerable 
differences appear in determining the matrix- valued Moyal bracket and in the classical limit. To consider this question 
one should find the matrix-valued Weyl symbol from the product of two operators in a standard way. 

Let A{p, q) and B{p, q) be matrix-valued Weyl symbols of two operators A and B respectively. Let us introduce 
operators Ci, (72,(7 as follows: 



Ci =AB 
C2 = BA 
C=[A,B]^Ci 



(11) 



C2 



Then, using 



and consideration like that used in [h5| , one can obtain the matrix- valued Weyl symbol of the operator 



Ci{p,q) = A{p,q)exp 



— {dqdp- dpdq) 



B{p,q). 



(12) 



For the matrix-valued Weyl symbol of the operator C'2 one can find a similar expression. Therefore, the matrix-valued 
Moyal bracket can be determined in the following form: 



{A{p,q),B{p,q)}M^^C{p,q)^ 



^(A(p,g)exp fCdqlfp-'dplfy) B{p, q) - B{p, q) exp f (d q'^ p - 1) plf q) A{p,q) 



(13) 



Unlike the usual case, matrix-valued Weyl symbols do not commutate, so it is not possible to represent such a Moyal 
bracket as a sinus from the operator of the Poisson bracket. This fact results in the classical limit in some peculiarities. 
In the general case, one can express the classical limit of the matrix- valued Moyal bracket via the matrix- valued Poisson 
bracket pl[ and the commutator of the matrix-valued Weyl symbols: 



lim{i,B}^, = 

^i^O iTl 



^{{A,B}p-{B,A}p) 



(14) 



If the matrix-valued Weyl symbols (for example position and Hamiltonian) commutate, the matrix-valued Moyal 
bracket is in fact the usual one and coincides with the Poisson bracket in the classical limit. If they do not commutate 
with each other, three cases are possible. Let [A, B] = O {fi^) when fi ^ Q. If j/ < 1, the classical limit does not exist 
for such a couple of the matrix-valued Weyl symbols, since the value of the matrix-valued Moyal bracket increases 
infinitely. If > 1, it becomes zero in the classical limit, li v = \, the classical limit exists. Here we can take as an 
example the Newton -Wigner position operator ||Tl|,|l^. Part of its matrix- valued Weyl symbol, not commutating with 
the Hamiltonian, is proportional to Planck's constant, thus in this case there exists a well-defined classical limit. 



III. MATRIX- VALUED WIGNER FUNCTION AND QUANTUM LIOUVILLE EQUATION 



Using (^ , one can find the average value of an arbitrary operator A expressed via its matrix- valued Weyl symbol: 

+ OC 

(15) 



A^ I A^^{P.q){^p\Wc.''{P.q)\'^")dpdq. 

^ Cl LI 



Q,/3,7=±l, 
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Then, introducing the matrix-valued Wigner function 



a=±l 

this expression can be written in a simpler form: 

+00 



A = tr J A{p, q)W{p, q)dpdq. (17) 

— C30 

Substitution of the operator of quasi-probability density expansion into ( |l6|) leads to the expression for the 
matrix-valued Wigner function in the coordinate representation. As a result we obtain a formula that coincides with 
the definition given in P0|: 



+00 



dQ. (18) 



Next, let us describe a method to obtain evolution equation for the matrix-valued Wigner function of a scalar 
charged particle. To do this one shall differentiate expression (^8|) with respect to time. The value of the wavefunction 
derivative is taken from the integral form of the Klein - Gordon equation (|^). After standard transformations, similar 
to the usual WWM formalism, we obtain the quantum Liouville equation 

dtW^{H,W}M, (19) 

where the matrix- valued Moyal bracket is determined by expression (^3|) . 

At first sight, this equation has a considerable defect in comparison to similar expressions in |l|,|l|; namely, it is not 
Lorentz invariant. For an interpretation of this fact, we have to take into account arguments from the measurement 
theory. The fact is that average values calculated in this approach coincide with ones in the usual (Schrodinger) 
representation of quantum mechanics that is Lorentz invariant. Nevertheless, unlike the stochastic formulation of 
quantum mechanics ||25[| , the scalar product is determined here with functions that are square integrable in a certain 
space-like hyper-surface, rather than over the whole space-time. Moreover, the value of scalar product does not 



depend on its choice |26|. One can relate this hyper-surface to the measurement device frame. In other words, the 
wavefunction collapse occurs in a frame in which equation ( p^ ) is written. The absence of Lorentz invariance is a 
consequence of the fact that the Weyl rule does not include time as an independent dynamical variable. 

In p7[ | the process of relativistic measurement is considered. In that work the point of view is expressed, that 
the wavefunction (and as result the Wigner function, we notice) has no objective value and does not covariantly 
transform when there are classical interventions. Here we note that in principle equation ( p^ can be written with 
four-dimensional Lorentz invariant symbols only, but to do this we have to incorporate in the theory a certain time-like 
unit vector in a way similar to Tomonaga-Schwinger approach to quantum field theory |^,^. It is the four- velocity of 
the frame where the wavefunction collapse occurs (the measuring device frame) relative to the second static observer 
(watching observer) . This frame has principally another sense to that of the preferred frame in [fzj-p^ . In that approach 
P,p^ collapse of the wavefunction has to take place in all frames in the whole space (because of the instant time 
hyper-surface being invariant). In our approach this process obeys the relativity of simultaneity. If it were possible 
to observe the wave function directly, such a hypothetical watching observer would see the wavefunction collapse as a 
certain moving front. As a result, at a certain instant, in a static frame (attached to the watching observer), a state 
with a part of the wavefunction before measurement on one hand and a part of the wavefunction after measurement 
on the other hand, is r ealized (Fig.|l]). But perhaps it is impossible to propose even gedanken experiment without use 
of a superluminar signal to interpret the results in favor of one of the approaches. 

Hence we find to be very important and of fundamental value the fact that quantum mechanics in the Wigner 
representation necessarily includes the four-velocity of the measuring device frame in explicit form. 
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Measuring device frame 
(moving observer) 




Static frame 
(watching observer) 
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FIG. 1. The process of the wavcfunctiori coUapsc in a moving frame from points of view of an observer in this frame and 
a static observer (watching observer). For the latter this process is a certain moving front. At a certain instant a state exists 
where at one side (right) there is the wavefunction before the measurement and on the other side (left) there is the wavefunction 
after the measurement. 



IV. TRANSFORMATION TO THE FESHBACH - VILLARS REPRESENTATION 



In the Feshbach - Villars representation, Hamiltonian matrix has a diagonal form and indices in the charge space 
correspond to the particle and antiparticle fl^ . Therefore, it is convenient to distinguish between solutions with 
different charge signs and to give them explicit physical sense. Furthermore, in this approach an influence of the odd 
part of the position operator on values of some physical variables is more evident. 

Operator Aq, ''transforms to the Feshbach - Villars representation by the following formula: 



\FV p 



where the transformation matrix has the form: 



-1 s 



7,<5=±1 



(20) 



[{E{p) + m^)Sj + {E{p) - m^)T,^P] 



(21) 



Here and below 



E{p) 



(22) 



is the energy of a relativistic free particle. 

Next, we apply (|2^) to the Weyl rule in the form (^). As result, we obtain the expression that brings into 
correspondence the matrix-valued Weyl symbol to the operator in the Feshbach- Villars representation: 



aFV 13 _ 



1 



p 

P+2 



E 

7,5=±1 



U^^iP + ^Wip. q)U-\\p - ^) exp 



/ P 
dPdpdq VP - 1^ 



(23) 
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In general, it is rather difficult to interpret this expression since there is a complicated dependence on the integration 
variable under the integral sign. Here we do not go beyond the simplest case. 
Let the matrix-valued Weyl symbol be proportional to the identity matrix: 



Aj{p,q) = A{p,q)dJ 



(24) 



In principle, one can say that such symbols do not depend on the charge variable, so that the class of dynamical 
variables, which corresponds to those, we denominate here, to be more brief, as a class of charge-invariant variables. 
Most of the dynamical variables that we consider in relativistic (non-quantum) mechanics belong to this class. The 
reason for this is the absence of a dependence on the charge variable in classical mechanics. Hence, the question about 
the classical limit is for such variables especially interesting. 

The Weyl rule for charge-invariant variables in the Feshbach - Villars representation has the form: 



A 



FV f) 



1 



(27r?i)° 



P\ P P 

A{p, q)Rj{p+ —,p~ ^)exp 



/ P 
dPdpdq ( P - Y 



(25) 



Unlike and non-rclativistic case there is a matrix-valued variable here: 

Rj{Pl,P2) = £{Pl,P2)Sa'^ + X{Pl,P2)Tla'^- 

It contains even and odd parts and is expressed via the energy of a free particle (|2 

e(pi,p2) = — f 

^^P^^P^) 2^E(p,)E{p2) 



(26) 



(27) 



Consequences from (^5|) are expressions for even [A\ and odd {^1} parts of the operator of a charge-invariant observable 
in terms of its Weyl symbol: 



A 



FV 



1 



(27ra)° 



A{p, q)e{p - 



-^,P- ^jOa exp 



-Pq 



I P 
dPdpdq VP - 1^ 



(28) 



+ OC 



/3 _ 



p\ p p 



/ P 

dPdpdq VP - 1^ 



(29) 



Matrix elements with eigenvectors of the momentum of the operator of a charge-invariant variable and its even and 
odd parts can be written in the following form: 



1 



+ 00 



{p^\A^'' J \P2) ^ j^—^Rc''{pi,P2) / A(-(pi+p2),9)exp 



,1 



-{pi ~P2)q 



dq, 



(30) 



{Pi 



A 



FV 



\P2 



+ 00 



{2TTh)' 



je{pi,P2)6a'^ 



A{^{pi +P2),q)exp 



-(Pi -P2)q 



dq, 



(31) 



+ 00 



(pil {i^^} J \P2) = j^x{pi,P2)riJ J A{ 



hiPi +l'2),9)cxp 



't;{pi ~P2)q 



dq. 



(32) 



Then, as was achieved in Sec.||, one can obtain a formula that reconstructs the Weyl symbol from the operator, with 
even and odd parts, in the Feshbach Villars representation: 
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A{p,q)6j= J2 

7=±1_ 



I P 



exp 



I 



dP, 



(33) 



_ 



-1/ ^ P^ 



^2 



~,FV 



(34) 



4-CXD 



7=±1_ 



P P, 



P + 



P 



P 



exp 



-Pq 



dP. 



(35) 



Comparing ( pl| ) and ( p2[ ) we conclude that matrix elements (integral kernels) of even and odd parts of the operator 
of an arbitrary charge-invariant variable are uniquely related to each other due to the Weyl rule: 



A 



FV 



\P2 



(36) 



A consequence of this expression is the fact that independent of position the odd part of an operator is zero. If one 
uses as A, for example, the scalar potential of an electric field, expression (^5|) establishes a quantitative relationship 
between effects of motion of a particle in an electric field and its interaction with a polarizable vacuum (trembling 
motion, Zitterbewegung) ]l2[ |. 

Consider now time derivative peculiarities of charge-invariant variables. In the coordinate representation, the 
matrix-valued Weyl symbol of such an operator has the form: 



dtA^{A,H} 



M- 



The Weyl rule (23) (in the Feshbach - Villars representation) for this can be written as follows: 



(37) 



dtA 



FV 



1 



or, in another form: 



^ \u{p + |){i(p, q), H{p)}mU-\p - |) exp 



-Pq 



I P 
dPdpdqip- — 



+ 00 



P + f ) {A{p,q)Mp' + f - ^)}^=''exp 



-Pq 



dPdpdq (p — 



(38) 



(39) 



where we have introduced a new matrix-valued variable of three arguments: 

E'ip) 



Gj {pi,P,P2) 



2^E{p,)E{p2) 

If A{p^ q) on q linearly depends, formula ( |39| ) can be presented in a particularly simple form: 

+00 



« aFV h f I X 9A{p,q) dE{p) o /3n , / I 

dtA'^^'a' = J \P) Q;^{T3a'^ +lT2a'^)dp{p\. 



(40) 



(41) 



The time derivative of an operator is a composition of its even and odd parts, where the odd part has a classical limit 
because of the absence of interference terms. As an example one can take the Newton - Wigner position operator 
that is even part of the usual position operator. 

In the general case this formula contains in higher orders of h non-standard terms, since G{pi,p,P2) distinguishes 
from E{p). However, in the classical limit they vanish and this expression takes the form of (4l[) 

Hence, in the Feshbach - Villars representation, where we distinguish between solutions wit 
the odd part of the position results not only in the emergence of the odd part of operators, but leads to some 
peculiarities in their even parts as well. In particular, it appears in the specifics of the semi-classical limit. 



1 different charge signs. 
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V. WIGNER FUNCTION AND QUANTUM LIOUVILLE EQUATION FOR CHARGE-INVARIANT 

VARIABLES 



It is easy to see from (|2^) that it is possible to introduce the usual Wigncr function for the charge invariant variables 
in such a way that their average values are determined by the formula: 



A{p, q)W{p, q)dpdq. 



(42) 



To show this, we expand the operator of quasi-probability density in the Feshbach - Villars representation through 
the eigenvectors of momentum: 



W'^Jip^q) 



1 



+ OC' 



(27r?i)" 



^\ nr P p. 
P+^ ) RJ {p+ ^,p- -^)exp 



(43) 



The Wigner function is determined as the average value of this operator over an arbitrary state that contains in the 
general case components with both charge signs: 



wip,q)= {Mw^''c.^ip,g)\r) 

a,l3=±l 



(44) 



There are four components in this expression. Two of them are the average value of the even part of the operator of 
quasi-probability density, and the other two are the average values of the odd part. Let us introduce the symbols: 



W^^{p,q) = {MW^''Jip,q)\r) 



(45) 



It should be noted here and below that the object Wa^{p, q) is not the matrix- valued Wigner function in the sense of 
M and SecJml of this work. 



Substituting ( |43| ) into ([45|), we obtain for the Wigner function components following expressions 

+ OC' 



1 



P P P P 



-Pq 



dP, 



(46) 



Wo.-"ip,q) 



P P P _ P 

X{P+ j,P- ^Wa{p+ ^)ip °'{p-j)exp 



dp. 



(47) 



Even components of the Wigner function (46) correspond to a charge definite state. The value of odd components 
(^) for such a state is zero. The expression (|4^ ) differs from analogous one for a non-relativistic Wigncr function 
and relativistic one determined using the Newton - Wigner position operator by the function e(pi,p2)under the 
integral sign (see (27)). This function has a specific feature: its expansion on pi, p2 does not contain square terms. 
This means that in the non-relativistic limit expression ( ^ ) coincides with the usual determination of the Wigner 
function. 

We obtain the evolution equations for every component separately. The general principle here is the same as in 
Sec.IIL However, instead of equation (||) we shall use the integral form of Klein -Gordon equation in the Feshbach - 
Villars representation [p3]: 



indt-^a.{p) = Eip)T3^c.-^f,ip). 



(48) 



The following equations can be obtained in a standard way, through differentiating Wigner function components with 
respect to time: 



dtWo.'^ip, q, t) = a'^E[p) sinl-^'a'p^ JW^,"(p, q, t), 



(49) 
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dtW^-'^ip, q, t) = ta^Eip) cosi^t pt.jWa^-'^ip, q, t). 



(50) 



Nevertheless, the Wigner function components are not independent, i.e. a specific constraint is imposed on solutions 
of the system (p9|),( p0[). To find this one should take the Fourier transform and make the standard change of variables 
for every component (|46|),(^). As a result, we obtain the following expressions: 



^*+{pi)-$+{p2)^e-\puP2) J W++{^{pi+p2),q)exp 



^{pi -P2)q 



dq, 



(51) 



i-OO 

\E'l(pi)*"(p2) = e"^(pi,P2) J VK_-(i(pi +P2),g)exp 



P2)q 



dq, 



(52) 



*;(P1)*-(P2) 



X ^{Pl,P2) 



w+ 



P2),q)e^p 



P2)q 



dq, 



(53) 



+ 00 

«'l(pi)*+(P2) = X"'(Pl,P2) J W- + {^{pi+p2),q)exp 



-{pi -P2)q 



dq. 



(54) 



Now we divide ( p2D by ( |53D and ( |54| ) by (51), and equate the resulting expressions with each other due to the equality 
of their left-hand sides. This gives us the constraint we are looking for: 



+ 00 



(S(pi) - E{p2)f J W++{^ipi +P2),qi)W^-i^{pi +P2),92)exp [^(gi + 92)(pi - P2)] ^91^92 = 

-°°+oo ■ (55) 

= {E{pi) + Eip2)y J W+-{^{pi +p2),gi)VK_ + (i(pi +P2),g2)exp [j-{qi + q2){pi - P2)] dqidq2 



Equation ( pOD explicitly contains the imaginary unit, so the actual question is whether the Wigner function is real. 
To test this we consider the complex-conjugate expressions for (|4^),(^7[). After some easy transformations we obtain 
the following identities: 



W\''{p,q) = W^"{p,q), 



(56) 



W^*a""(p,9) = W^-a"(p,g). (57) 

These mean that even components of the Wigner function are real. Odd components are complex conjugate to each 
other so that their sum is real as well. 

It is essential that the equation for even components of the Wigner function coincides with the analogous expression 
obtained in p3[ | for the formalism where the Newton - Wigner position operator is used. Hence, the dynamics of quasi- 
distribution functions for systems of particles with charges of the same sign (charge definite states) is identical in both 
cases. 



VI. STATISTICAL PROPERTIES OF THE WIGNER FUNCTION FOR CHARGE-INVARIANT 

VARIABLES 

Constraint on the initial conditions of the Wigner function is the general peculiarity of the approach described here 
because equations are identical in both cases (for charge definite states) . In this Section we show how some theorems 
and properties differ from their analogues in the usual WWM formalism Jl6[ and in an approach where the Newton - 
Wigner position operator is used p^. 
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First of all one should note the property of normality. Even part of the Wigner function ( fig ) is normalized in the 
whole phase space, and the integral of the odd part is zero. 

Consider now the compatibility of the Wigner function ( ^6| ) with distributions in the coordinate and momentum 
spaces for a single-charge state. For this purpose we integrate (|4^ ) by coordinate. As a result, we obtain the 
distribution in momentum space: 



(58) 



This function always has a definite sign, so that it can be interpreted as the probability density. One can obtain a 
more non-trivial result for distribution in coordinate space. The result of integrating of (46) by momentum can be 
given in the following form: 



or, which is the same, 



1 



+CC 



*a(Pi)exp 



e(pi,P2)*"(p2)exp 



dpidp2 



(59) 



(60) 



where 'i>a{q) is the wavefunction in the representation of the Newton - Wigner coordinate It is obvious, that 



quasi-distribution, (g9|) and 
probability interpretation [BO 



(60D, is not sign-definite. This property is typical for bosons and causes difhculties in 
a%. Nevertheless, formal use of such quasi-probability makes it possible to calculate 
average values of the variables dependent on coordinate. 

Average values of variables only depending on momentum do not differ from similar ones in usual approach. Let 
us find the peculiarities of higher moments of coordinate in our case. For this purpose we integrate (^9|) with q" by 
coordinate and use the identity: 



g"gi(p2-'Pi)g 



{—ihdp^Y 6xp 



^(P2 -pi)q 



(61) 



After some obvious transformations the result for the n-th moment of the coordinate can be written as follows: 



'I 

*"(p)£(p,p') dp. 

J p'=p 



(62) 



The first moment (average coordinate) has a value similar to one in the Newton - Wigner coordinate approach. 
Differences manifest themselves in higher moments. As an example, we apply this expression to second moment of 
the coordinate: 



+ 00 



+ 00 



(9') = / Kip) 



ih d , 



Kip) 



hc^p 



2E^p) 



^°'{p)dp. 



(63) 



As well as the usual part, this formula also contains an additional term that causes the peculiarities related to 
determination of the position operator. For strongly localized states, it results in formal violation of the uncertainty 
relation. 

Similar to the case of the usual WWM formalism ||l^ , we prove two criteria that make it possible to select Wigner 
functions for pure and mixed states out of the whole set of functions of the variables (p, q) . 



Criterion 1 Criterion of pure state. 

For the Junctions Wa'^{p,q) and Wa~^'^{p,q) to he even and odd components of the Wigner function for charge 
invariant variables, it is necessary and sufficient that equalities (|5^, ([5^,([5^ hold true, and the following conditions 
are satisfied: 



+00 



dpidp2 



■In 



Wa°'{^{pi +P2),g)exp 



-{Pi 



P2)q 



dq 



c^PlP2 



E{pi)E{p2){E{pi) + Eip2))^' 



(64) 
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+ 00 



^{pi ~P2)q 



dq 



c'^PiP2 



E{pi)E{p2){Eipi) - Eip2W 



(65) 



Necessity. Following ||l6|] we start from an obvious identity: 

92 



ln«':(pi)*'3(p2) = 0. 



dpidp2 

Applying it to formulas (51) - (54) one obtains the equalities 



+00 



-^In / M^a"(kpi+P2),'?)exp 
opidp2 J 2 



^(pi -P2)q 



dq = 



dpidp'. 



■lne(pi,p2), 



(66) 



(67) 



+00 



^In/W^«-"(^(P1+P2),g)exp 



^{pi -P2)q 



dq 



dpidp: 



■lnx{pi,P2)- 



(68) 



Now, substituting the explicit form of e{pi,p2) and x(piiP2)fron i (|27|), w e obtain (p^ , (p^) . 

Sufhciency. Let components of the function satisfy conditions (64),(|65|) or, equivalently, (|6^),(p8|). So one can write 
following conditions: 



+00 



e ^(Pi,P2) / W^a"(-(pi +P2),g)exp 



-{pi -P2)q 



dq =(pi{pi,a)ip2{p2,a.), 



(69) 



+ 00 



X ^{pi,P2) I Wc "(-(pi +p2),g)exp 



^(pi -P2)q 



dq ^ip3{pi,a)ip4,ip2,a), 



(70) 



where ipi{pk,a) are certain functions. Let us show that they can be chosen in such a way as to be consistent with 
each other and present a wavefunction for a scalar charged particle. From ([5^),(|69|) and (|57|) , ([70|) it follows that 



ipl{pi,a)ip2{p2,a) = ipi{p2,a)ip2ipi,a), 



(71) 



Then one obtains 



ipl{pi,a)ipl{p2, a) = (P4(P2, -a)ip3{pi, -a). 



f2{pi,a) if^{p2,a) 



= const = A, 



const = B. 



tp4{pi,a) tpl{p2,-a) 
Hence, one can say that our system is described by the following wavefunctions: 



Kip) 



VA ' 



(72) 

(73) 
(74) 

(75) 
(76) 



^-"ip) = VBMP,a), 



(77) 
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(78) 



Now one has to prove that wavefunctions presented in such a way (with one and two tildes) are consistent with each 



other. For this purpose, a combination of expressions ( |75| ) - (78) and ( |69| ) and (|7C|) is substituted into condition (|5 
As a result, we obtain 

I;(pi)f-(P2)^-(P1)^+(P2) = *;(Pl)*~(P2)*-(Pl)*+(P2)- (79) 

This leads to the following equality: 

nipi)^-ipi) ^-(P2)^+(P2) , ^ 

= = = const = 1. (80) 

*^^(pi)*l(pi) *-(p2)*+(P2) 

If this condition is fulfilled, one can accept as the wavefunction 

^^{p) = l„(p). (81) 

Then condition ( [70| ) is satisfied and, taking account of the fact that components of ^'^(p) differ from components of 
^aip) by a constant phase, condition (|69| ) is satisfied as well. Hence, conditions of this criterion allow us to introduce 
wavefunction (pT|). Then, similar to how it is done in ||l^, one can show that it satisfies the Klein - Gordon equation 
in the Feshbach - Villars representation. Which was to be proved. 

This criterion differs from a similar one in the usual WWM formalism [|l^ and in the approach that uses the 
Newton - Wigner position operator by difference of the right-hand parts of conditions (|6^),(|65|) from zero. The 
obvious consequence from this is the fact that the one-particle state, where Wigner function would be a joint Gauss 
distribution by coordinate and momentum, is impossible in the approach described here. 

Consider another property of Wigner function, true for both pure and mixed states. For this purpose, we introduce 
the formula for the square of the module of the scalar product of state |\E') and state |<i>): 

\{^\^)\^ ^ {27rhy E J W^^^{p,q)e-Hp+^^d„p+^^d,)W^„''{p,q)dpdq+ 

Tc^'-"° ■ (82) 

+ i27rhY E / W^o.~"ip,q)x-^ip+^,'dg,p+^^tg)W^-^''ip,q)dpdq 

a— ±1 — oo 

It is easy to test this by substituting expressions (ji^ ) and (^) for the Wigner function components into it. Then, let 
W{p, q) presents a mixed state consisting of orthogonal states described by Wn{p, q)- 

WJip, q)^Y. ''nWnJip. g), (83) 



Taking into account that the sum of 5n(the probabihty of system to be in state n) squares is less than unity (a 
consequence from ( ^4| ) ) one obtains for an arbitrary state 



E / W^"{p,q)e-^{p+^^d„p+^^d,)W^''{p,q)dpdq+ 

"=^'-°°+oo ■ (85) 

+ E J W^'"{p,q)x-HP+Trd,,P+0,)W-a''{p,q)dpdq<j^ 

Moreover, for a pure state this inequality turns into an equality. Expression (|8^) can be used as necessary and 
sufficient condition for both pure and mixed states. 

This condition is written for a state that is the superposition of states with different charges signs, which is more 
typical for the many-particle case. It takes a simpler form for a state where only one charge sign is realized: 

J Wo."{p,q)e-^{p+ ^t„p+ ^'^,)Wo."{p,q)dpdq < j^^. (86) 
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Though this inequahty does not contain interference terms, it differs from its non-relativistic analogue due to function 
(^) being present. Hence, a dynamical variable that contains square and higher moments of coordinate can have 
peculiarities in many-particle systems that are described by statistical physics. 

The type of influence of the odd part of the position operator on dynamical variables can be illustrated with the 
one-particle state that is the Gauss distribution in momentum space with arbitrary square of dispersion Ap^: 



: exp 



27rAp2 



4Ap2 



(87) 



where we have chosen the natural units of measure, h — c = m = 1. The dependence of position dispersion on 
momentum dispersion is shown in Fig.^j. Its peculiarity is the formal violation of uncertainty relation. Under very 
strong localization, states even with negative square of dispersion are possible, which was shown in |30 3^. This fact 
is a property of scalar charged particles, and localization peculiarities for fermions rather different. 




Ap 



0,06" 

Ax 

0.045 - 



I ' I ' I ' I ' I ' I 

6.6 6.8 7.0 7.2 7.4 7.6 

Ap 

FIG. 2. Dependence of the momentum dispersion on the position dispersion for the "Gauss" state. Dispersion of the 
coordinate is given in Ac — units; dispersion of the momentum is given in mc units. For comparison the curve AX = 
for usual uncertainty relation is shown. 



VII. CONCLUSIONS 



The usual (not Lorentz invariant) Weyl rule makes it possible to introduce the Wigner function that is not Lorentz 
invariant, but all average values calculated with it coincide with ones calculated with Lorentz invariant wavefunction. 
This results in the fact that quantum mechanics in the Wigner formulation contains with necessity a measuring device 
frame. In principle, we can write the evolution equations using only four-dimensional Lorentz invariant symbols, but 
it is necessary to introduce a certain time-like vector for it. It is the four- velocity of the frame in which wavefunction 
collapse occurs, relative to the second (static) observer. It should be noted, that this approach differs from [p|jic|] 
where there is the preferred frame. In the last cases such a frame has a global sense and its introduction is related 
to attempts of correct tachyons description and, as a consequence, to a possible explanation of the instant quantum 
correlations (in a relativistic case) from the position of de Broglie - Bohm quantum mechanics. 

Phase space for a scalar charged particle is not only limited by three couples of the momenta and coordinates . 
The charge part of it exists as well. However, in the approach presented here we leave the operator nature of such 
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variables without modification. As a result, the matrix- valued Wigner function is the density matrix in charge space 
with standard rules of average values calculation as well. 

If we limit our consideration only to such elements of dynamical algebra that do not depend on variables of the 
charge space, it is possible to introduce the usual Wigner function. This object differs from the Wigner function for 
a non-relativistic particle and from the Wigner function in the Newton - Wigner position operator approach as well. 
First of all it should be noted that it contains four components, corresponding to particles, antiparticles and their 
interference with each other. Moreover, even for the one-particle case, when only one component exists, the definition 
of the Wigner function differs, as the result of odd part of the position operator being present, from the usual one. 

This results in non-standard behaviour of some physical variables, even in the absence of conditions when particles 
creation is possible. However, it is not observed for all physical variables. For example, energy and number of particles 
(that are usually considered in statistical physics) do not show such peculiarities. Hence, one can expect such effects on 
the quadratic and higher moments of coordinate. As an example, one can take the dispersion that can be interpreted 
as the real (physical) size of system that is not limited by external borders. 

One can separate two groups of effects that result from this approach: ones related to interference between particles 
and antiparticles, and ones that take place in systems with same charge signs. Effects of the first group result from 
presence of the odd part of the Wigner function, the second ones are due to the specific function e{pi,p2) being 
present in the even component of the Wigner function. At the one-particle level it manifests itself, for example, in 
violation of the uncertainty relation. Perhaps, such effects can exist in many-particle systems as well. 

Even and odd components of Wigner function in the system of quantum Liouville equations are not mixed up 
together. This results from absence of particles creation from vacuum in the system. For example, in an electric 
field the vacuum is not stable |3|], and this can be interpreted as a consequence of the odd part of the position 
operator being present as well. It manifests itself in mixing up of different components in the quantum Liouville 
equation. There is another situation for the instant magnetic field: particles are not created and the Wigner function 
components are not mixed up |2^ . Nevertheless, even components of position and momentum of such a system satisfy 
not the usual commutation relations but the deformed Heisenberg - Weyl algebra ones |Q. These facts can mean 
that in external electromagnetic fields the odd part of the position operator reveals itself to be especially strong. 
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